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ABSTRACT 

Aims. We present X-ray and l2 CO(J=l-0) observations of the very-high-energy (VHE) y-ray source 
HESS J 1813—178 with the aim of understanding the origin of the y-ray emission. 

Methods. High-angular resolution X-ray studies of the VHE y-ray emission region are performed using 
18.6 ks of XMM-Newton data, taken on HESS J1813-178 in October 2005. Using this dataset we are able to 
undertake spectral and morphological studies of the X-ray emission from this object with greater precision 
than previous studies. NANTEN 12 CO(J= 1-0) data are used to search for correlations of the y-ray emission 
with molecular clouds which could act as target material for y-ray production in a hadronic scenario. 
Results. The NANTEN 12 CO(J=1-0) observations show a giant molecular cloud of mass 2.5 x 10 5 M© at a 
distance of 4 kpc in the vicinity of HESS J 18 13— 1 78. Even though there is no direct positional coincidence, 
this giant cloud might have influenced the evolution of the y-ray source and its surroundings. The X-ray 
data show a highly absorbed (n H ~ 1. x 10 23 cirT 2 ) non-thermal X-ray emitting object coincident with the 
previously known ASCA source AX J 18 13— 178 showing a compact core and an extended tail towards the 
north-east, located in the centre of the radio shell-type Supernova remnant (SNR) G12.82-0.2. This central 
object shows morphological and spectral resemblance to a Pulsar Wind Nebula (PWN) and we therefore 
consider that this object is very likely to be a composite SNR. We discuss the scenario in which the y-rays 
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originate in the shell of the SNR and the one in which they originate in the central object in terms of a 
time-dependent one-zone leptonic model. We demonstrate, that in order to connect the core X-ray emission 
to the VHE y-ray emission electrons have to be accelerated to energies of at least 1 PeV. We conclude that 
if indeed the X-rays are connected to the VHE y-rays HESS J18 13-178 has to be a Galactic Pevatron. 



Key words. ISM: Supernova remnants, plerions - ISM: individual objects: HESS J1813-178, G12.82- 
0.02, AX J 1 8 1 3—178 - X-rays: observations - gamma-rays: observations 



1. Introduction 

During a survey of the Galactic plane in very high-energy (VHE) y-rays the atmo- 
spheric Cher enkov telescope system H.E.S.S . found more than a dozen new >100 GeV y- 



ray sources (lAharonian et al 



2005a, 



2006al) . Even though some of the new y-ray sources 



can be identified with counterparts a t other wavebands, for example the microquasar 



LS 5 039 (lAharonian et al 



part dFunk 



2005b 



2006bl) . most of them lack a compelling positional counter- 



2006b - A firm counterpart identification for a y-source at other wavebands requires 



not only a spatial association, but also a viable y-ray emission mechanism and a consistent multi- 
frequency picture. X-ray observations with instruments such as XMM-Newton and Chandra com- 
bine high-angular-resolution and high-sensitivity observations with the possibility to trace rela- 
tivistic electrons via their synchrotron emission. These instruments are thus ideally suited to 
the identification of VHE y-ray objects. This article reports on multi-frequency observations of 
HESS J 1 8 1 3—178 with the XMM-Newton X-ray satellite. Dense molecular regions that could act 
as target material for hadronic interactions have been studied with the NANTEN CO-telescope 
in the direction of HESS J 1 8 1 3—178 to test for an origin of the y-ray emission in the decay of 
neutral pions produced in proton-proton collisions. 

HESS J 1 8 1 3—178 was in fact the first source discovered during the survey. It was subse- 
quently re-observed to increase the statistical significance and the signal in 9.7 hours of live- 
time comprised 340 y-ray events at a significance of ~ 14 standard deviations. While many 
of the new sources found in the survey showed a rather large source extension, HESS J 1 8 1 3— 
178 has a compact nature and is only slightly extended (Gaussian width cr = 2.2' + 0.4') with 
respect to the point-spread function (PSF) of the H.E.S.S. instrument. The flux of the source 



is (14.2 + l.lstat) X 10 photons cm s above 200 GeV and the energy spectrum follows 
a power-law with photon index 2.09 + 0.08. HESS J 1 8 1 3—178 is located in the Galactic plane 
at 18 h 13 m 38 s , -17 d 50'33". HESS J1813-178 was subsequently also detected by th e MAGIC 



Cherenkov telescope confirming the basic properties of the source dAIbert et al 



2006) 



2005al) . HESS J1813- 



At the time of the first publication of this source ( lAharonian et al. 
178 was still considered as unidentified. However, subsequently to the H.E.S.S. publication, the 
source was reported to show compelling positional coincidence with several sources in other 
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energy bands dBrogan et al 



2005 



Helfand et al 



2005 



Ubertini et al 



2005). HESSJ1813-178 



coinci des with the previou sly unpublished archival ASCA source AX J18 13-178 in the 2-10 keV 



band dBrogan et al 



2005). This unresolved ASCA source with a rather hard photon index of 
1.8 + 0.4 and a flux of 7 x 10~ 12 erg crrT 2 s~' (2 -10 keV) was one of the relatively bright sources 



2001). AXJ1813-178 exhibits a 



detected in the ASCA Galactic Plane survey (ISugizaki et al. 
highly absorbed (n H = (1.1 ± 0.2) x 10 23 cm- 2 ) non-thermal spectrum without indications of 
emission lines. The value of nH is significantly larger than the total column density through 
the Galaxy in this direction, suggesting that the source is embedded in a dense environment. 
Positionally in coincidence is the unresolved hard X-ray source IGR J18 135-1751 detected i n 



2005) 



the INTEGRAL Galactic plane survey in the 20-100 keV energy band dUbertini et al 
The source connects to the energy spectrum of the ASCA source, showing however a softer 
energy spectrum, suggesting a spectral break in the X-ray spectrum around 10 keV. Recently the 
source has also been observed with the Swift X-Ray Telescope (XRT) that found an X-ray source 
coincident with the ASCA source. The XRT data confirm th e hard spectral ind ex and the high 



column density towards the source present in the ASCA data (iLandi et alj|2006l) 



The region surrounding HESS J 1 8 1 3—178 was also covered by observations in the 20 cm and 
90 cm band by the VLA. These data sho wed a faint non-thermal radio source (G12. 82-0. 02) 



visible at the position of HESS J18 13-178 dBrogan et alj|2005 



White et al.ll2005l) . G12.82-0.02 



lies in the projected vicinit y (distance ~ 15') of the W 33 region. This region is known to contain 
ultra-compact HII regions ( Church wel J 1990 ) and contains methanol, hydroxyl, and water masers 
and other tracers of recent star formation. While H.E.S.S. ASCA, and INTEGRAL lack the 
spatial resolution to resolve the object, the angular resolution of the VLA allows a shell-like radio 
structure of diameter 2.5'to be resolved. The flux at 20 cm was reported as 0.65+0. 1 Jy, the flux at 
90 cm was given as 1 .2 ± 0.08 Jy, yielding a non-thermal radio spectrum of index a = -0.42+0.0 3 



(j v oc v"). Spitzer Space telescope data from the GLIMPSE survey at 8 fim dBenjamin et al 
2003b show no signs of dust emission in positional coincidence with G12.82-0.02, a finding that 
supports the non-thermal natu re of the source. Th e radio spectrum, morphology and the lack 



of coinciding IR-Emission led IBrogan et al 



(12005b to the conclusion that the radio structure is 



a previously unknown young shell-type Supernova remnant (SNR). In the Parkes multi-beam 
pulsar survey (PMBS) no pulsar close to G12.82-0.02 has been found thereb y determining the 



1.4 G Hz flux density upper limit at a rather constraining level of 0.2 mJy dManchester et al 



200 lh . Currently the only constraints on the distance of the object can be derived from the strong 
absorption found in the X-ray data that suggests that AX J 1 8 1 3—178 is located behind W 33 at a 
distance larger than 4 kpc. 

Assuming an association between HESS J 1 8 1 3—178 and the radio source G12.82-0.02, the 
question remains whether the thus far unresolved X-ray and y-ray emission originates from 
the s hell of the SNR, or rather from a Pulsar Wind Nebula (PWN) embedded within the 
shell. IBrogan et al 



(2005) concluded from the fact that the radio spectrum nearly directly ex- 
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trapolates into the ASCA X-ray spectrum, that both the radio shell and the ASCA source should 
indeed be the same source and therefore that all the emission was connected to the radio shell. 

This article reports on high-angular resolution X-ray observations performed with XMM- 
Newton with the goal of pinning down the origin of the high-energy X-ray and y-ray emis- 
sion. Another option for the y-ray emission is the interaction o f accelerated hadro ns with dense 



molecular material. The bright ultra-compact HII region W 33 (IChurchwellll 1 9901) is the nearest 
star forming region, located ~ 10' away from HESS J 1 8 1 3—178. Dense gas in the W33 region 
could act as target material for the VHE y-ray generation in hadronic interactions, therefore a 
study of the 12 CO(J=1-0) dis tribution performed during the NANTEN survey of the Galactic 
plane ( Mizuno & Fukuill2004 ) towards HESS J1813-178 is included. 



2. XMM-Newton Observations of the region 

HESS J1813-178 was observed with XMM-Newton on the 14fh of October 2005 for 18.6 ksec. 
The instrument cameras (EPIC MOS1, MOS2, PN) were operated in full-frame mode with a 
medium filter to screen out optical and U V light. The data were processed with the XMM-Newton 
Science Analysis Software (SAS) version 7.0 as well as with the Extended Source Analysis 
Software package (XMM-ESAS) version 1.0. Standard data reduction and calibration proce- 
dures were applied to the data and the total data set amounts to 13.6 ksec after the observations 
have been cleaned of flaring background caused by soft protons. Figure Q] shows the resulting 
combined count map of the EMOS1 and EMOS2 detectors for three different energy bands (red: 
0.5 keV to 2 keV, green: 2 keV to 4.5 keV, blue: 4.5 keV to 10 keV) smoothed with a Gaussian 
kernel of 3 pixels. No shell-like emission as seen in the VLA data but rather a compact X-ray 
source with an extended tail towards the north-west is visible at the centre of the EMOS-Cameras. 
This source is in positional coincidence with the previously reported ASCA source AX J 1 8 1 3— 
178. There is still a contamination of soft photons in the dataset. This had been seen before in 
the ASCA data and can be attributed to the bright, close-by (~ 0.7°) low mass X-ray binary 
GX 13+1. To investigate the influence of the soft photon contamination in the background re- 
gions for the spectral studies, three separate methods were applied to estimate the background as 
discussed later in the text. The results for the different methods show good agreement, enhancing 
the confidence that the background level in the spectral studies was estimated correctly. There is 
however still the possibility of soft photon contamination in the source spectrum, especially at 
energies below 2 keV, since the source appears strongly absorbed below this energy. 

For the source position fitting the co ntamination was avoided by using the source detec- 
tion algorithm (emldetect) as described in ISnowden. Collier & Kunta d2004l) for energies above 



4.5 keV and for energies above 7.5 keV. Seven sources are detected in the band above 4.5 keV 
which are not in the part of the field of view affected by the stray light (see Table [T), whereas 
above 7.5 keV only a single X-ray source is detected. This source is extended and coincident with 
AX J1813-178 at a best-fit position 18 h 13 m 35.16 s , -17 d 49'50.0", with a statistical error on the 
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Table 1. X-ray sources other than AXJ1813-178 detected in the field of view surrounding 
HESS J 1 8 1 3—178 using the detection algorithm emldetect. The second column gives the name 
recommended by the XMM-Newton SOC and the IAU for serendipitous XMM-Newton source 
detections. Column 3 and 4 give J2000 coordinates with a typical error on the position of 1". 
Column 5 gives the number of counts in a circle of radius 30"above 2 keV. 



Id 


XMMUJ 


RA2000 


Dec 2 ooo 


Counts 






(h:m:s) 


(d:':") 




1 


181259.2-175527 


18:12:59.23 


-17:55:26.9 


60 


2 


181314.1-175344 


18:13:14.11 


-17:53:44.1 


238 


3 


181318.0-175232 


18:13:17.57 


-17:52:31.9 


90 


4 


181323.7-175041 


18:13:23.68 


-17:50:41.2 


138 


5 


181348.3-175301 


18:13:48.35 


-17:53:00.6 


68 


6 


181333.3-175857 


18:13:33.27 


-17:58:56.7 


157 




AXJ1813-178 


18:13:35.16 


-17:49:50.0 


2192 



position of 2". This position is at a distance of ~ 1' from the best fit position of HESS J18 13-178, 
located well within the extension of 2.2' of HESS J 1 8 1 3—178. Even though an a priori match of 
the best fit X-ray and y-ray positions is far from being certain as e.g. shown in the PWN candi- 
date HESS J 1825-1 37 (in whic h the H.E.S.S. position is significantly shifted from the peak of 



the X-ray emission as shown in 



Aharonian et al 



d2006d) ) this X-ray position is compatible with 



the H.E.S.S. position, given a statistical error of 21" adding to a systematic error in the pointing 
accuracy of the H.E.S.S. instrument of the similar size. The source position fitting tool emldetect 
also determines AX J 1 8 1 3—178 to be incompatible with a point-source. Using a Gaussian model 
the extension was determined to be 8.5 + 0.2 image pixels, corresponding to an projected width 
on the sky of 21" + 0.5". As apparent in the slice through the source, shown in Figure |2j the 
Gaussian profile might not be the correct representation of the source profile, which seems to 
have a rather long tail towards the north-east. The source position determined using the XMM- 
Newt on data is compat ible with the position of the only X-ray source found in the Swift XRT 



data (ILandi et al 



20061) at 18 h 13 m 34.9 s , -17 d 49'53.2"with an error of 3.5". The distance of the 



best fit position presented here to the XRT position is 4.9". The other X-ray sources within the 
field of view listed in Table[TJare no known X-ray sources. 

For the spectral analysis XSPEC (version 12.2.1) was used and three different background 
methods as describe above were applied to estimate the stray-light contamination on the spec- 
trum of AX J18 13-178. The background was estimated a) from the south-western part of the field 
of view, where the stray-light contamination is apparently lower b) from a ring around the source 
w ith inner radius 100"and outer rad ius 200"and c) from blank field observations as described 
in lSnowden. Collier & Kunta d2004l) . All three methods yield consistent results for the spectral 
fit, in the following the background derived from a ring around the source within the same field 
of view will be used. Different extraction radii of size 50", 75", and 100" were used to determine 
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Table 2. Spectral properties of AXJ1813-178 as determined by the spectral fitting routine. 
EMOS1, EMOS2, and EPN data were simultaneously fit with a single absorbed power-law. The 
fit parameters are the absorption density nn, the photon index F, and the norm at 1 keV, which is 
expressed by the flux between 2 keV and 10 keV. Different spectral shapes such as a black-body 
spectrum were incompatible with the data (just to give an example: the^ 2 /d.o.f. for a blackbody 
spectrum for the extraction region of size 75"amounts to 5700. The errors given here correspond 
to the 90% confidence levels. 



Radius 
(arcsec) 


Counts 

EMOS1 EMOS2 EPN 


n H T F 2 -iokeV ^ 2 /d.o.f. 
10 22 citT 3 10~ 12 erg cnT 2 s" 1 


50 
75 
100 


2004 1805 3930 
2957 2769 6194 
3794 3592 8397 


12.0^;°. 1.81!^ 5.2 3177/3001 
10.6^ 1-72+H;! 2 6.8 3233/3001 
10.6!°j 1.881°;°| 69 3228/3001 



the spectrum of the extended emission. All EMOS1, EMOS2, and EPN data were fit simultane- 
ously. In general the data are well described by a single power-law and the results of the spectral 
fitting are summarised in Table [2] Figure [3] shows the XMM-Newton spectrum for the medium 
size extraction radius of 75". As can be seen, no prominent line emission can be found and anal- 
ysis shows that the shape is incompatible with a black-body radiation spectrum, confirming the 
non-thermal nature of the X-ray emission. The results for the different extraction radii are well 
compatible with each other in terms of column density and spectral index. The column density 
determined is «h ~ 10 23 cm~ 2 , significantly hig her than the total Galactic column density in this 



region of the Galactic plane, 1.6 x 10 22 cm 1 (iDickey & Lockman 



19901) . The photon index is 



determined to be ~ 1.8, the flux between 2 and 10 keV for the 75" extraction radius amounts 
is ~ 7 x 10~ 12 erg cirT 2 s _1 . Both these values are well compatible with previous measurements 
reported fro m the analysis of ASC A data (photon index: 1.8 + 0.4, F2-10 keV = 7.0 x 10~ 12 erg 



cm 2 s 1 ) by Brogan et al 



(120051) . 

The EPN data within a narrow region (5") surrounding AX J 1 8 1 3—178 was searched for pe- 
riodic emission from any possible X-ray pulsar. The timing resolution of the EPN detector in 
full-frame mode is 73.2 ms. No significant periodicity was found. Also no evidence for variabil- 
ity was detected in this dataset. Future X-ray timing observations with XMM-Newton or Chandra 
and future Multibeam radio observations might shed further light on the existence possible un- 
derlying pulsar. 



3. NANTEN-Observations towards HESS J1 81 3-1 78 

The 12 CO(J=1-0) observations were performed with the 4-meter mm/sub-mm telescope, 
NANTEN, at Las Campanas Observatory of Carnegie Institutions of Washington, Chile. The 
NANTEN 12 CO(J=1-0) survey of the Galactic plane shows that the velocity of the 12 CO emis- 
sion towards HESS J18 13-178 ranges from vLSR (velocity local standard of rest, i.e. the average 
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velocity of stars in the solar neighbourhood) ~ km/s to ~ 60 km/s. While the CO emission at 
vLSR less than 10 km/s is likely due to local clouds, the emission at velocities in excess of 
vLSR ~ 20 km/s is likely from clouds beyond 2 kpc. The best known object in the region sur- 
rounding HESS J 1 8 1 3—178 is W33, a complex of HII regions, a typical example of a massive 
star forming region. The integrated intensity distribution at vLSR — 30^40 km/s shows an en- 
hancement in close projected vicinity to H ESS J1813— 178 as sh own in Fig. [4] To convert the 
velocity range into a distance the model by iBrand & Blita d 19931) has been used and the rather 



large uncertainty in kinematic distance for directions close to the Galactic centre has to be taken 
into account. The W 33 complex has been studied fairly extensively in the past and its distance 
has bee n estimated to be ~ 4 kpc (kinematic distance derived from molecular line observations 



see e.g. dReifenstein et al 



4 kp 
1970: 



Goldsmith & Mao 



1983 



Mitchell et al 



1990)). The most mas- 



sive molecular cloud in the direction of HESS J 1 8 1 3—178 is in the velocity range 30-40 km/s as 
shown in @] The size of this cloud is ~ 70 pc x40 pc in Galactic longitude and latitude (at an 
assumed distance of 4 kpc) and its total mass is estimated to be 2.5 x 10 5 solar masses from the 
12 CO(J=1-0) velocity integrated intensity, equiv alent to 10 21 cm' 2 in molecular column density, 



assuming an X-factor of 2 x 1 20 cm 2 / (Kkms 1 ) (iBertsch et al 



19931) (where X = JV(H 2 )//(CO), 



is the conversion factor between H2 column density and integrated CO intensity). This mass is 
typical for a giant molecular cloud. HESS J18 13-178 is clearly located outside of the dense cloud 
on a scale of ~20 pc and therefore the cloud and the VHE y-ray emission are probably unrelated. 
The region surrounding HESS J 1 8 1 3—178 is a fairly complicated ensemble of young massive 
stars as indicated by several c ompact HII regions. Fi gure[5]shows Spitzer Space telescope data of 



the GLIMPSE survey at 8 fim ( iBeniamin et al 



20031) . As previously mentioned, no dust emission 



at the position of HESS J1813-178canbe found (red contours denote the 1, 2, and 3 cr error con- 
tours on the centroid of the VHE y-ray emission), however, the peak in the CO-data seems to be 
correlated to the 8 //m data, supporting a star-formatting origin of the giant molecular cloud. The 
fact that still a reasonable amount of dust in the region that has not yet been blown away by the 
hot stars might suggests that these stars are still rather young. The total integrated column den- 
sity determined for a beam size of 2.6'centred on Galactic longitude 12.8, latitude 0.0 (the closest 
position to HESS J1813-178 in the NANTEN s urvey with a grid size of 4') is ~ 9 x 10 22 cirT 2 , 



somewhat higher than the value determined by 



Dickey & Lockmanl d 19901) (1.6 x 10 22 cm 2 ), 



but fully consistent with the high absorption density determined from the XMM-Newton data. 
Integrating the column density only to 4 kpc, the proposed location of AX J18 13-178 results in a 
valu e of ~ (0.5 - 1.) x 10 22 cm~ 2 . Comparing the NANTEN values presented here with the values 



from 



Dickey & Lockmanl (119901) the differences are due to both the difference in material probed 



(atomic versus molecular hydrogen) as well as the angular resolution of the two surveys. 
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4. Interpretation of the X-ray emission, connection to other wavebands 

Even though the NANTEN CO-Emission is probably unrelated to HESS J1813— 178, it is inter- 
esting to calculate an upper limit on the cosmic ray density in the dense molecular cloud given 
the non-detection of a VHE y-ray signal. To calculate the ratio k of the cosmic ray density in the 
cloud pcr.gmc to the local cosmic ray density pcR.iocai the following relation is used: 

F(>\ TeV) X Z) 2 

k = PCR.GMC / PCR.local = — — ~~~ 7Z ~ " 5 TTT ( 1 ) 

2.85 x 10~ 13 x M 5 (kpc- cirr 2 /s) 
(here F(>i TeV) is the y-ray flux above 1 TeV in units of cur 2 s" 1 , D is t h e dist ance in kpc, 



Aharonianl (119911) . The factor 



and M5 is the mass of the cloud in units of 10 5 solar masses) 
2.85 x 10~ 13 kpc 2 crrT 2 /s has been derived assuming a CR spectrum with a photon index of 
2.6. Using a 2cr H.E.S.S. upper limit of 0.6% of the Crab flux above 1 TeV (taken from the 
10 hours exposure time and the H.E.S.S. sensitivity of 1% Crab flux in a 25 hours observation), 
i.e. 1.1 x 10~ 13 cirT2 s , a distance of D — 4 kpc, and a cloud density of M5 = 2.5, an upper 
limit onp CR; GMC ~ 3 xp CR ,iocai can be derived. 

Figure |2] shows a comparison between the X-ray and Radio emission in the region surround- 
ing HESS J 1 8 1 3—178. As can be clearly seen, the radio emission shows a shell-like structure 
whereas the X-ray emission has a compact core with an extended em ission towards the north- 
east, a typical morphology for a PWN (see iGaensler & Sland (120061) for a recent review). The 



apparent anti-correlation with the radio shell also suggest a confinement of the X-ray emission 
within the shell, especially since the tail of the X-ray source extends to the NE, where a break 
in the radio shell is present. The non-detection of a pulsar in reasonably deep radio observa- 
tions may be due to beaming effects. In the X-ray data presented here the thermal emission 
from the possible neutron star might be buried underneath the strong non-thermal emission and 
additionally strongly absorbed by the high column density which strongly affects the detection 
of X-rays below ~ 1.5 keV. However, to finally confirm this scenario, the pulsar within this 
nebula would need to be found, either in deep radio or X-ray observations. Alternative strong ev- 
idence that AX J 1 8 1 3—178 is a PWN would be the detection of spectral softening away from the 
cor e, a signature of ele ctron cooling which is observed in many PWN, see for example G21.5- 



0.9 



Slaneetal 



(120001) . The XMM-Newton data do not allow such an effect to be resolved in 
AX J 1 8 1 3—178. A search for a steepening in the data with respect to the distance of the best fit 
centroid yielded no conclusive results within statistical errors. Nevertheless, it seems very likely 
that AX J 1 8 1 3—178 is indeed a PWN, as the positional coincidence of a hard spectrum X-ray 
source within a radio shell of G12.82-0.02 is otherwise very unlikely. Furthermore, a chance 
positional coincidence of the VHE y-ray source HESS J 1 8 1 3—178 with this composite radio/X- 
ray object also seems unlikely, an in the following discussed we will assume that the emission 
seen in these three wavebands originates in a single new composite SNR, similar in its properties 



to e.g. GO.9+0.1 (lAharonian et al 



|2005cJ). However, the situation in the case discussed here is 
somewhat different, since the VHE y-ray source shows extended emission (possibly owing to the 
smaller distance of HESS J1813-178 (~ 4 kpc) in comparison to GO.9+0.1 (~ 8 kpc)). 
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Both the shell of the SNR as well as the newly found PWN candidate AXJ1813-178 are 
viable y-ray emitting objects. To estimate whether the shell of the SNR could be responsible for 
the y-ray emission the following prescription has been followed: a background (determined in 
a ring of 2—\' around the SNR centre) has been subtracted from the central 2' radio emission, 
the resulting radio map has been smoothed with a Gaussian kernel of width r smoot h = 1.2', slices 
have been fit through the resulting smoothed map in RA- and Dec-Direction with Gaussians to 
determine cr x and <r y of the smoothed emission region. Finally the "Gaussian equivalent width" 
of the SNR has been defined as cr SNR = (cr 2 + cr 2 - r 2 mooth ) 1/2 , yielding cr SNR ss 1.8', compatible 
with the measured rms extension 2.16' + 0.36' of the VHE y-ray emission region. The SNR shell 
seems therefore to be a plausible candidate for the origin of the y-ray emission. 

As a result, we have to conclude that even with high angular resolution XMM-Newton X-ray 
observations in which the central X-ray emitting object AXJ1813-178 was resolved, no final 
distinction can be drawn between the scenario in which the y-rays originate in the shell or in the 
core of G12.82-0.02. The X-ray data described here and previously described radio observations 
indicate that G 12.82-0.02 is a composite supernova remnant with a bright X-ray core and a radio 
shell. The size of the gamma-ray source measured by H.E.S.S. appears to be consistent with an 
origin of high energy emission in the SNR shell, but a common origin of the X-ray and y-ray 
emission in a central PWN cannot be excluded as a larger spatial extent of the y-ray source with 
resp ect to the X-ray source c ould occur in such cases (see for example the case of HESS J 1825- 



137 ( lAharonian et al 



2006cJ)), due to the energy dependent cooling of electrons. Two distinct 



scenarios for the origin of the TeV emission must therefore me considered: 1) as a counterpart to 
the X-ray emitting core, and 2) as a counterpart to the radio em i ssion of th e shell. 
Scenario 1) has previously been discussed by 



Ubertini et al 



(120051) and 



Helfand et al 



(2005) 



Here the situation is revisited in the light of the new XMM-Newton data. Figure [6] (top) shows 
the spectral energy distribution of G 12.82-0.02. The EGRET upper limit has been derived from 
the first 5 years of the EGRET mission yielding a flux upper limit of 2.7 x 10~" erg cirT 2 
s _1 above 100 MeV. Where angular resolution is sufficient (i.e. in the radio and < 15 keV X- 
ray bands) the core and shell of the remnant are shown separately. Two synchrotron/Inverse 
Compton model curves are shown for a population of relativistic electrons in the core. The model 
was chosen to be time-dependent with constant injection over the lifetime of the source. The 
key model parameters are the slope of the injection spectrum of electrons a the minimum and 
maximum energies of the electrons E m \ n and E m . dx , the magnetic field in the source B, and the 
target radiation field (CMBR, optical or dust photons). In both c ases shown an age of 1000 



years is adopted, in-line with the estimates of 



Brogan et al 



(120051) for the age of the remnant: 



300-3000 years. The combined XMM- Afewton/INTEGRAL spectrum indicates that E max must be 
rather high: ~ 10 15 eV. If inverse Compton emission takes place in the Thompson regime then 
equal keV synchrotron X-ray and VHE y-ray IC spectral indices are expected away from any 
cut-off in the electron spectrum. The softer spectrum measured by H.E.S.S. therefore suggests 
that Klein-Nishina (KN) effects may be important in this source. The well-detected X-ray peak 
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determines roughly a convolution of the magnetic field with the square of the maximum electron 
energy B x E mRX 2 , whereas the ratio of the total synchrotron to IC emission determines the ratio of 
B 2 to the radiation density, i.e. the magnetic field B if assuming normal radiation fields. Therefore 
the electron index a is the only free parameter in the modelling when assuming normal radiation 
fields. A rather soft electron spectrum is needed to match the y-ray data, whereas the X-ray data 
suggests a harder spectrum. The dashed lines in Figure [6] (top) show a model with a = 2.4, 
B = 4.2 fiG, E m i„ = 25 GeV, £ max = 1.5 Pe V and a radiation field with near and far infra-red 
components each of energy density 1 eV cm 3 in addition to the CMBR. Such a scenario is in 
marginal agreement with the H.E.S.S. and X-ray data, but a low energy break or cut-off in the 
electron spectrum is required to avoid over-producing radio emission in the SNR core. This cutoff 
could be explained by the termination shock of a PWN, in which the electrons in the lab-system 
have gained a minimum energy through bulk motion. This minimum energy can be as high as 
~ 100 GeV up to 1 TeV, for a typical PWN Lorentz factor between 10 5 and 10 6 . The spectral 
break can be avoided by invoking a strong contribution of scattering by IR/optical radiation fields, 
where IC scattering transits into the KN regime and is thus cut off at higher energies, resulting 
in an apparent steepening of the VHE y-ray spectra compared to the X-ray spectra. The solid 
curves in Figure [6] (top) show a model with a — 2, B — 7.5 fiG, E mm = 1 MeV, E max = 1.5 PeV 
and a radiation field with a very strong NIR component, with energy density 1000 eV cm" 3 . This 
radiation field exceeds b y a factor of » 1000 the nominal NIR radiation field at 4 kpc from the 
Galactic Centre (see e.g. 



Porter et al 



(2006)) and is thus somewhat unrealistic. The (possibly) 



nearby star formi ng region W 33 may contribute to the radiation density in the vicinity of the 



Helfand et al 



(12005b have estimated the radiation density in G 12.82-0.02 to 



source. However, 

be 3-4 eV cm -3 , including the contribution of W 33. Sce narios with intermediate val ues provide 



Hint on & Aharonianl (120061) for details 



equal acceptable agreement with the available data. See 
of the calculation methods used here. 

A leptonic model, in which the core X-ray and VHE y-ray emission are associated yields an 
unavoidable £ max > 1 PeV, suggesting that HESS J 1 8 1 3—178 is a highly effective accelerator - a 
Galactic Pevatron. In case of a soft electron injection spectrum E max could well be higher, since 
the XMM-Afewto/i-INTEGRAL spectral break could plausibly be a cooling break, rather than 
represent the end of the electron spectrum, for a somewhat greater pulsar age. If the X-ray and 
VHE y-ray emission are indeed associated (i.e. originate from the same electron population), 
the PeV maximum energy for the accelerated electrons implies the emission almost certainly 
has plerionic origin. PWN are certainly capable of accelerating electrons to PeV energies as 
e.g. shown by EGRET measurements of the Crab Nebula, showing a synchrotron component 
extending to MeV energies (shell type SNRs in the framework of diffusive shock acceleration 
scheme fall short of 1 PeV for electron acceleration). The total energy injected in the electrons 
is 6 x 10 46 ergs for the high-IR case and even 6 x 10 47 ergs for t h e soft electron spectrum case. 



For a source of age 1000 years (as suggested by iBrogan et al.l (120051) ) these values result in 
an electron luminosity of L e ~ 2 x 10 36 - 2 x 10 37 ergs/s. This value can be compared to the 
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comp arable composite SNR GO. 9+0.1 which has a derived electron luminosity of L e ~ 7 x 10 37 



erg/s dHinton & Aharonianl2006l) 



A distinct alternative scenario is an origin of the y-ray emission in the SNR shell. Two young 
shell-type SNRs (with estimated ages comparable to that of G 12.82-0.02) are well established 
VHE y-ray sources, namely RXJ17 13.7-3946 and RX J0852.0-4622 (known as Vela Junior). 
Placed at the greater distance of 4 kpc, the TeV luminosity of HESS J18 13-178 is comparable to 
that of these two SNRs. Furthermore, the measured size of the TeV emitting region is consistent 
with an origin in the radio shell. The SNR shell must therefore by considered seriously as an 
alternative source of the TeV emission. Figure [6] (bottom) shows two models for a y-ray origin 
in the SNR shell: a) as inverse Compton emission from the electron population responsible for 
the radio emission (solid line), and b) as the product of the decay of neutral pions produced in 
hadronic interactions of acc elerated protons in or near the SNR shell (dashed line), calculated 



using the parametrisation of 



Kelner et al 



d2006l) . In either case, the central X-ray emission has 
to attributed to a different mechanism, presumably a pulsar. The electron model has parameters: 
a — 2.1, B — 5 yuG, E m \ n = 1 MeV, E max = 30 TeV and a nominal 4 kpc radiation field (note that 
larger values of E mm up to 1 GeV are still compatible with the radio emission). The small value 
of £max = 30 TeV is required to avoid producing significant X-ray synchrotron emission in the 
shell, for values much larger than this, 2 keV emission from the shell should have been seen in 
this XMM-Newton observation. The proton model has a = 2.1, E m \ n = 1 GeV, £max = 100 TeV. 
To match the TeV flux level in the p-p interaction scenario the product of the total energy in 
protons and the ambient density must be: (£ , p /10 50 erg)(«/l cm -3 ) » 1. For likely densities 
greater than 1 cm -3 , the required acceleration efficiency is < 10% for a typical SNR explosion 
energy of 10 51 erg. In case of a very large local density, Bremsstrahlung effects start to become 
important and might increase the VHE y-ray flux in comparison to the synchrotron X-ray flux. 
However, as we consider most of the absorption in the X-ray spectrum to be foreground density, 
the density within the source is likely insufficient for Bremsstrahlung to be the dominant process. 



5. Summary and Conclusion 

Detailed XMM-Newton X-ray and NANTEN 12 CO(J=1-0) studies have been performed towards 
the y-ray source HESS J 1 8 1 3—178. The NANTEN data show a giant molecular cloud in the 
vicinity of HESS J 1 8 1 3—178 that might have played an important role in the evolution of the 
y-ray source. The X-ray data show a resolved non-thermal object in the centre of the shell of the 
radio SNR, most likely representing synchrotron emission from a PWN. The confirmation of this 
picture would require the detection of a central pulsar, either in X-ray or in radio observations, 
or indirectly through the detection of spectral cooling as seen in several other PWN systems. 
Nevertheless, G 12.82-0.2 is now tentatively established as a composite SNR through its distinct 
morphology of a central extended non-thermal object within a radio shell. For the first time, an 
object of this type has been detected first in VHE y-rays and then identified with high-angular 
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resolution radio and X-ray data. This detection shows, that y-ray observations are well suited 
to identify SNR, that are otherwise very hard to detect due to obscuration. No distinction is so 
far possible between a scenario in which the y-rays are emitted from the shell of the SNR, and 
one in which they are emitted from the central PWN. An upcoming deep Suzaku exposure on 
HESS J 1 8 1 3—178 will shed more light on the situation in the hard X-ray band and the future 
GLAST satellite will provide important constraints in the MeV-GeV band. If the central X-ray 
source and the VHE y-ray source are indeed connected our modelling of the SED emission 
suggests that HESS J1813— 178 is a Galactic Pevatron. 
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Fig. 1. Top: A composite false-colour image of the EM0S1 and EMOS2 count maps, slightly 
smoothed with a Gaussian kernel of width 3 image pixels. Different colours correspond to 
0.5 keV-2 keV (red), 2 keV-4.5 keV (green), and 4.5 keV-10 keV (blue). AXJ1813-178 ap- 
pears as a bright, hard-spectrum source with a tail towards the north-west. The north-western 
edge of the field of view shows some soft stray-light contamination, caused by the low mass X- 
ray binary GX 13+1. Also shown are several other X-ray sources that were found with emldetect 
at energies above 4.5 keV. These are typically fainter and have a softe r photon index. The white 



Brogan et al 



(120051). The 



contours denote the VLA 20 cm radio emission as already reported by|] 
shell-like radio structure surrounding AX J 1 8 1 3—178 is clearly visible. Bottom: Radial profile of 
the central source AX J 1 8 1 3—178 in comparison to the PSF for this data set. 
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Fig. 2. Comparison of radio and X-ray data of HESS Jl 8 13-178. Top: XMM-Newton counts 
map above 4.5 keV of the region surrounding HESS J18 13-178 (colour contours) smoothed with 
a Gaussian kernel of width 0.1'. The extended tail towards the north-east is visible in this figure 



Overl aid is the 20 cm shell-like emission (white contours) as detected by the VLA (iBrogan et al 



20051) . The difference in the images between the X-ray and the radio wavebands is apparent in 



this figure. Also shown are t he positional contours ( 1, 2, 3 <x error) of the best fit position of 



HESSJ1813-178 as given in 



Aharonian et al 



d2006al) (dashed circles) and the extension (solid 
green), covering both the radio and X-ray emitting region completely. Bottom: Slice through the 
emission in radio and X-rays as plotted on the left hand side. The box in which the slices were 
determined is also given in the left panel (white box). The X-ray slice shows the compact core 
with the slice towards the north-east, whereas the radio slice shows the shell-like structure of the 
emission. 
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Fig. 3. XMM-Newton energy spectrum for all three detectors (EMOS 1 : black, EMOS2: red, EPN: 
green) for an extraction radius of 75", rebinned to yield a minimum significance of 5cr in each 
bin, summing a maximum of 100 bins into a single one bin. The straight lines show the combined 
fit to these data (folded with the instrumental response which is different for the three detectors). 
The lower panel shows the residuals of the data to the fit. As can be seen, the power-law fit 
provides a reasonable description of the data. Also apparent is the absence of prominent line 
emission that lends support to the non-thermal nature of the emission. 
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Fig. 4. NANTEN 12 CO(J=1-0) data towards the region surrounding HESS J18 13-178 for three 
different velocity ranges. Please note that this figure is in Galactic coordinates. The colour scale 
shows the CO emission in units of K km s _1 , the red contours the H.E.S.S. VHE y-ray excess 
emission (smoothed with a Gaussian of 0.05° width and p lotted in equidistant excess contours), 



whereas the white contours show the 20 cm VLA contours (IBrogan et al 



20051) . The dashed cyan 



box indicates the size of the area shown in Fig. [2] The shell-type radio emission at the centre of 
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Fig. 5. Spitzer Space telescope data of the GLIMPSE survey at 8 pm (Beni amin et al. 
in Galactic coordinates. The red contours show again the H.E.S.S. VHE y-ray excess emission 
(smoothed with a Gaussian of 0.05° width and pl otted in equidistant excess contours), whereas 



2005b . The green contours show 



the white contours show the 20 cm VLA contours (Br ogan et al. 
the NANTEN CO-Emission contours at the distance of 4 kpc, showing the position of the GMC. 
The dashed cyan box indicates the size of the area sho wn in Fig. [2] No dust emission is coincident 
with the VHE y-ray source as reported already by iBrogan et al.l (I2005I) . supporting the non- 
thermal nature of the source. The Spitzer and NANTEN data are clearly correlated on this scale. 
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Fig. 6. Spectral energy distribution for HESS J1813-178. The XMM-Newton data 2 keV and 
10 keV are shown for an extraction radius of 75" and are corrected for the absorption. A system- 
atic error band of 0.2 on the spectral in dex and 20% on the fl ux level has been added. For the 
INTEGRAL points the data published in 



Ubertini et al 



(2005) have been recently reanalysed to 



determine error bars on the flux points. The flux points shown here are consistent with the pub- 
lished po ints. The H.E.S.S. spectr um has been rebinned in comparison to the energy spectrum 
d2006al) . Top: Scenario in which the VHE y-ray emission originates 



shown in 



Aharonian et al 



in the core of the SNR (i.e. the PWN). Leptonic one-zone model (lAharonian & Atovan 



1999) 



for two scenarios in which the VHE y-ray emitting and the X-ray emitting electrons belong to 
the same population. The model parameters are specified in the text. Bottom: Scenario in which 
the VHE y-ray emission originates in the radio shell of the SNR. Leptonic (solid) and hadronic 
(dashed) VHE y-ray emission scenarios in which the y-ray emission originates in the shell of the 
SNR, i.e. the radio and VHE y-ray emitting electrons are connected. 
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